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ABSTRACT 

We describe a model for the thermal and dynamical equilibrium of starless 
cores that includes the radiative transfer of the gas and dust and simple CO 
chemistry. The model shows that the structure and behavior of the cores is 
significantly different depending on whether the central density is either above 
or below about 10^ cm^^. This density is significant as the critical density for 
gas cooling by gas-dust collisions and also as the critical density for dynamical 
stability, given the typical properties of the starless cores. The starless cores thus 
divide into two classes that we refer to as thermally super-critical and thermally 
sub-critical. This two-class distinction allows an improved interpretation of the 
different observational data of starless cores within a single model. 

Subject headings: CO 

1. Introduction 

Starless cores are dense regions (riHj ~ 10^ to 10^ cm^'^) in dark clouds with linear 
scales of tenths of pc, and total masses of a few solar masses. The starless cores contain 
no infrared sources above the sensitivity level of the IRAS satellite (about 0.1 at the 
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starless cores (lAlves. Lada. &: Ladall200ll : iTafalla et al.ll2004l ) show density profiles that ap- 
proximat ely match th ose of pressure-confined, hydrostatic spheres (Bonnor-Ebert or BE 
spheres ; iBonnorl Il956l ) . 



Despite an overall similarity in structure, recent observations with increasingly better 
angular resolution and sensitivity show that the starless cores are not all the same. For exam- 
ple, observations of L1517B indicate i sothermal gas and s how spectral line profiles indicativ e 
of little motion or possibly expansion (iTafalla et al.ll2004j : iKeto et al.ll2004l : ISohn et al.ll2007l ). 
In contr ast, observations o f L1544 indicate a significant variation in temperature from center 



to ed ge (ICrapsi et a 



tion (jWilliams et al 



■20071) and show spectral line profiles cons i stent with an ove rall contrac- 
I999I : ICaseUi et aljl2002l : keto et al.l[20oi : ISohn et~aDl2007h . The cores 



B68 and TMC-1 provide another example in contrast. Obser vations of CO in B 68 indicate 
an excitation temperature rising toward the edge of the core (iBergin et al. I2OO6I), and show 



spectral line profiles consistent with internal oscillations (so und waves) ( 



In comparison, observations of CO isotopologues in TMC-1 C (jSchnee et al. 



ada et all 120031 ). 



20071 ) indicate an 



excitation temperature decreasing toward the edge and show spectral line profiles consistent 
with contraction. 



In an earlier paper (IKeto fc Fieldll2005l ) we suggested that the observed differences were 
consistent with the different temperature and density structures of starless cores that have 
central densities greater or less than about 10^ cm~^. The significance of this particular 
density is twofold. 

First, this density divides the starless cores into those that are dynamically stable (those 
with lower densities) and those that must be contracting or collapsing (those with higher 
densities). The dynamical stability, of course, also depends on the total mass and the internal 
energy of the core, but for starless cores with masses of 1 - 10 Mq, temperatures of around 
10 K, and subsonic internal velocities, the critical density for dynamical stability falls within 
a range around 10^ cm~^. 

Second, at densities below a few 10^ cm~'^ the molecular gas in starless cores cools 
primarily through molecular line radiation, but at higher densities the gas cools primarily 
through coUisional coupling with the dust. This difference in cooling results in different 
temperature structures in cores with central densities below or above 10^ cm~^. If the density 
in a core is everywhere below 10^ cm~^ so that the gas and dust are nowhere coUisionally 
coupled, then the entire core cools primarily by molecular line radiation. As a result, the 
core is approximately isothermal. If the core has a central density above 10^ cm~'^, the gas 
in the core center cools by collisional coupling with the dust whereas the gas elsewhere cools 
by molecular line radiation. This difference in cooling results in a variation in temperature 
across the core. Because the dust, heated by external starlight, is cooler in the core center. 
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the gas there, coupled to the dust, is also cooler. 



In this paper we continue our investigation of the structures of the two classes of starless 
cores. We improve our previous model of thermal equilibrium by accounting for the heating 
of the gas by hot electrons photoelectrically released from dust grains and by including 
variable molecular abundances based on a simple model for CO chemistry. Our improved 
model for starless cores allows a better interpretation of the observational data. 

For example, the inclusion of photoelectric heating provides a better description of the 
thermal structure near the core boundary where this effect is most sign ificant. Applied to 



recent observation s of the edges of the starless cores B68 and TMC-IC ( iBergin et al.l 12006 



Schnee et al.ll2007l ) the rising and falling CO excitation temperature seen in the observations 
of the two cores, respectively, are consistent with the different effect of photoelectric heating 
on cores with relatively higher and lower densities at their boundaries. 

The loss from the gas phase of CO and other molecular coolants has on ly a minor 
effect on the gas temperature and hence the density structure of starless cores ([Goldsmith 
2OOII ). Nonetheless, an understanding of the molecular abundances is quite important in 
the interpretation of spectral line observations, particularly the line profiles and strengths. 
We have used our simple model of CO chemistry in modelling spectral line profiles in an 
oscillating lower mass, lower density starless core (iBroderick et al.l 120071 ) and in modelling 
spectral line strengths in a collapsing, higher mass, higher density, starless core (E. Keto & 
P. Caselli, in prep.). 

Our improved models of starless cores are also useful in constraining some of the prop- 
erties of the interstellar medium (ISM) that relate to the thermal equilibrium of the star- 
less cores. In particular, we obtain temperatures closer to those indicated by observations 
with dust opacities higher than those of naked grains and m ore consistent with fluffy dust 
(lOssenkopf fc Henninglll994l : iKruegel fc Siebenmorgenlll994j ). and cosmic-ray ionization rate 
towards the lower end of the range of estimates, 1 x 10~^^ s~^ to 6 x 10~^^ s^^. We find that 
photoelectric heating improves the stability of cores by increasing the confining pressure with 
reduced overlying mass, but additional internal energy is still required if the larger cores are 
to be supported. 
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The two classes of starless cores 



2.1. Observational description 

The higher mass and higher density cores are unstable to collapse and the formation 
of protostars within a free-fall time, but being starless, these cores have not yet formed a 
protostar. Thus these cores may be thought of as "young brides" that are soon to give birth 
to stars. The density profiles within these cores are more centrally concentrated owing to the 
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consistent with inward motions and 



line widths that increase toward the centers of the cores consistent with inward acceleration 
of the infall. Spectral lines from some cores sh ow two componen ts separated by more than the 



sound speed indicating supersonic velo cities (ISohn et al. 



2007). These denser cores are not 



isothermal . Observations of b oth dust ( Ward-Thompson. Andre, fc Kirkll2002l ; iPagani et al. 
2OO31 . 120041 : Tschnee et al.ll2007[ ) and gas JCrapsi et al.ll2007l : IPagani et aLlboOTn indicate cooler 



temperatures in their centers, the latter consistent with the change from molecular line cool- 
ing to more efficient dust cooling as the density increases toward the interior. 

In contrast, lower mass and lower density cores are stable against gravitational con- 
traction, and therefore in the absence of changes in their environment, will remain starless. 
These barren cores may be thought of as "old maids". These cores have flatter density 
profiles consistent with the relatively greater importance of an external pressure in confining 
the cores. Their spectral line profiles may be simple Gaussians indicating little i nternal mo- 



tion o r show complex shapes consistent with internal oscillation s (sound waves) (iLada et al. 



2003: Redman et al. 2006 



Aguti et al.ll2007l : iMaret et al.ll2007l). Multi- wave length infrared 



observations also indicate internal density perturbations ( Steinacker et al. 2005al lbl). The 
linewidths are generally more constant across the core and indicate subsonic internal veloci- 
ties. These cores are more nearly isothermal (~ 10 K) consistent with molecular line cooling 
throughout. 

Some recent papers have also drawn attention to differences within the population of 
starless cores and use terms such as "pre-stellar" to indicate a tend ency toward collapse and 
star formation versus "starless" or "failed" to suggest the opposite (I Andre. Basu &: Inutsuka 
20081 ). Sometimes the distinction is drawn as "gravitationally bound" and "gravitationally 
unbound" suggesting a ratio of greater or less than unity of gravitational energy to internal 
energy including thermal energy and possibly microturbulence and magnetic energy. These 
terms suggest a division of the starless cores into two classes distinguished by their gravi- 
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tational stability and future evolution. iKirk et al.l (120051 ) divide the starless cores into two 
categories called "bright" and "intermediate". These refer to the brightness of the dust 
emission which is related to the mass and density of the cores. In this paper we suggest that 
the distinction is motivated by other observational properties as well, and that the various 
properties of the two classes may be derived from a single theoretical model for the starless 
cores that we describe in the following several sections. 



2.2. Theoretical description: Thermally Super-Critical and Thermally 

Sub-Critical 



We refer to the two classes of cores as "thermally super-critical" and "thermally sub- 
critical". We choose these terms to emphasize that the distinction is based on the physics 
of the cores rather then their origin or evolution. The term "thermal" in this description is 
appropriate because it refers to both their thermal-radiative equilibrium and also to their 
thermally-supported dynamical equilibrium. Whether the cores are supported primarily by 
thermal energy or by magnetic energy is still a matter of debate. The evidence is gathering 
on the side of thermal sup port with non-thermal energy being an important but minor con- 



tribu t or to their stabil ity (IBarranco fc Goodman! Il998l : iGoodman et al.lll998l : iTafalla et al. 



2004 : iLada et al.ll2008l ). In other words, without any non-thermal energy most of the cores 
would be unstable. However, to achieve stability, the cores require only an amount of non- 
thermal energy that is a fraction of their thermal energy. The terms thermally super and 
sub-critical are similar in meaning to the terms magnetically super and sub-critical that are 
well established in the literature, but the "thermal" modifier suggests that thermal energy 
is here more important than magnetic energy. It is not necessary to prove this supposition 
to adopt the terminology, and in any case, in this study we model the cores as thermally 
supported. 

As typical parameters for our calculations, we adopt a central density of 10^ cm~^ and 
a total mass of 10 Mq to represent a thermally super-critical core, and 10^ cm~^ and Mq for 
a thermally sub-critical core. The we l l-stud ied cores L1544 and B68 are estima ted to have 
simila r cen tral densities. ICrapsi et al.l (120071 ) estimate 2 x 10^ cm"^ for L1544 and lKeto et al. 



( 2006 ) and Bergin et al. ( 20061 ) estimate 2.5 x 10^ cm ^ and 3.0 x 10^ cm ^ respectively for 



B68, all three estimates based on fitting isothermal BE spheres to the observational data. 
While L1544 and B68 are prototypes of cores in the two classes, here we are not attempting 
to model these two cores exactly. Our example parameters are chosen as round numbers that 
provide a factor of ten difference in density and mass, sufficient to bring out the differences 
in their structures. Of course, the real cores have combinations of mass and density that 
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form a continuum that includes these two examples. Other parameters for the cores are 
listed in table 1. 



3. CO chemistry 

Molecular abundances in starless cores are the result of complex interactions involving 
gas-phase chemistry, freeze-out and desorption, and photodissociation. In order to develop 
a simplified model we assume that the abundance of CO is decreased by two effects: 1) 
freeze-out onto dust grains in the high density center, and 2) photodissociation at the edge 
of the core. 



3.1. The photodissociation region 



Tielens fc HoUenbachl (Il985l ) suggest that at the edge of a molecular cloud where pho- 
todissociation is important, the dominant cycle for the formation and destruction of CO 
is, 

C+ CH2 CO ^ C ^ C+ (1) 

The time scale for the formation of CH2 by radiative association is longer than for the 
formation of CO from CH2 and O. Therefore we assume that the rate of formation of CO 
is given by that of CH2. This eliminates consideration of the complex chemistry of oxygen. 
The CO cycle may then be described by three rate equations for the creation and destruction 
of C"*", C, and CO, for example, 

(^C+ 

= -Rn{C^) + Rc{C^) (2) 

where Rn and Rc are the rates for the destruction and creation of C"^. Similar equations 
hold for C and CO. The conservation equation in the form, 

C+ + CO + C = 1 (3) 

indicates that the symbols stand for the non-dimensional relative abundance of the three 
species with respect to the total abundance of carbon. Assuming steady state we may solve 
for the relative abundance of CO as, 



r^r^ ^C0\ / CO C\ 



(4) 
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The ratios of the relative abundances, CO/C^ and C/C^, are equal to the ratios of the rates 
of destruction and creation of the respective species in the numerators. Using the data in 
Tielens &: HoUenbachl (Il985l ). these ratios are, 



and 



CO 

c+ 

c_ 

c+ 



1.4 X 10-"Goexp(-3.2Av 
2.1 X 10-iOGoexp(-2.6Av 



n(H2)6 X 10 



-16 



(5) 



(6) 



2.1 X 10-iOGoexp(-2.6Ay) 

Here, n(H2) is the number density of hydrogen, Gq is the interstellar radiation field in 
units of the Habing flux (Habing 1968), Ay is the mean visual extinction to the interstellar 
radiation field, and th e nurn erical values of the coefficients are taken from tables 5 and 12 of 
Tielens fc HoUenbachl (119851 ). The mean extinction at each point in the cloud is computed 
as the extinction from the point to the cloud surface averaged over all directions. 



(exp(-Ay)) 



1 

4:71 



exp{—Av)dQ 



(7) 



The Habing flux Gq = 1 corresponds to the average interstellar radiation field of the Galaxy. 
Figure [1] shows the variation of the three species, CO, C, and C^ as a function of Ay for 
a cloud with a constant density of 1000 cm~^ and Gq = 1. For this fiux level, the simple 
model provides for exponential rep lacement of CO and C "* " at a n Ay of about 1. The more 
complex photodissociation model of lTielens fc HoUenbachl (119851 ) shows that this exponential 
replacement is the dominant relationship, although in their figure 9b the replacement occurs 
at a higher Ay of 3 — 4 because of the higher radiative flux in their model. Go = 10^. 



3.2. Freeze-out 
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suggest that the molecular abundance in dense gas is an equilibrium between the rate of 
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Pavlyuchenkov et al.ll2003l : IShematovich et al.ll2003l : iLee. Bergin fc Evansll2004j : lAikawa et al. 
20051 ). This suggests that the steady state abundance of CO in dense gas may be computed 
as the ratio of the depletion time to the sum of the depletion and desorption times. 



The time scale for depletion onto dust may be estimated as (IRawlins et al.lll992l ). 

Ton = {SoRdgn{Y{^)aVT)-' (8) 

Here 5*0 is the sticking coefficient, with 5*0 = 1 meaning that the molecule sticks to the dust 
in each collision; R^g is the ratio of the number density of dust grains relative to molecular 
hydrogen; a is the mean cross-section of the dust grains; and Vt is the relative velocity 
between the grains and the gas. If the grains have a power law distribu tion of sizes with th e 
number of grains of each size scaling as the -3.5 power of their radii (IMathis et al.l 119771 ). 
then we can estimate their mean cross-section as, 



02 



n{a)da 



ax 



0.2 



n{a)a{a)da 



(9) 



ai 



where ai and 02 are the minimum and maximum grain sizes. If ai = 0.005 /xm and 02 = 0.3 
/im, then (a) = 3.4 x 10~^ fim^. Similarly, the ratio of the number densities of dust and gas 
may be estimated by computing the mean mass of a dust grain and assuming the standard 
gas to dust mass ratio of 100. If the density of the dust is 2 grams cm~^, then the ratio of 
number densities is Rdg = 4 x 10^^°. The relative velocity due to thermal motion is. 



Vt 



8kT 

TT/i 



'1/2 



(10) 



where T is the temperature and /z the molecular weight. 



There are several processes that generate heat on dust grains that can liberate frozen 
molecules to the gas phase. These include exothermic H2 formation on the grains, cosmic ray 
collis ions with dust grains, and photo-desorption of molecules off the grains (IWillacy fc Williams 
I993I ). The UV radiation for this la tter process de r ives fr om the ionization of the molecular 
gas by cosmic rays. According to [Roberts et al.l (120071 ). in starless cores these processes 
are all independent of density. The rate of H2 formation is independent of the molecular 
density because the density of HI in starless cores is independent of the molecular density. 
The rates of all these desorption processes are not well define d, nor is it known which is 
dominant. Therefore, for definiteness we adopt the rate given by lHasegawa fc HerbstI (119931 ) 
for desorption by cosmic rays but with the understanding that other processes might also be 
important. The time scale is. 



Toff = 3.3 X 10^ 



C 
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(11) 
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where the cosmic-ray ionization rate C = 3 x 10~^^ s~^, and the binding energy of CO onto 
ice is Eco = 1100 K (lOberg et aLll2005l . 120071 ). We choose this value for the binding ene rgy 
because it is between t he binding energie s of CO on water ice, 1180 K (IColhns et aLll2003a| ]bl). 



and on CO ice, 850 K fiFraser et aLll200 
of ~30 on ice mantles in dense clouds (jWhittet et al.ll2007l ). 



), assuming a typica l solid H2O/CO abundance ratio 



The relative abundance of CO in the gas phase in equilbrium between depletion and 
desorption is. 



CO 



gas 



Ton + Toff 



(12) 



where CO^a^ + CO grains = 1, and the symbols stand for the non-dimensional relative abun- 
dance of CO in the two phases with respect to the total CO abundance. The depletion time 
scales inversely with the collision rate while the desorption time scales inversely with the 
cosmic-ray ionization rate. Therefore the depletion of CO is dependent only on the gas den- 
sity because the cosmic-ray ionization rate is assumed to be the same throughout the cloud. 
Figure [2] shows the dependence of the steady state abundance of CO versus the molecular 
gas density. 

While we only illustrate the steady state abundance of CO in this paper, the model for 
depletion and desorption may also be used to describe the time dependent evolution of CO 
between the gas and grain phases. The change with time of the abundance in the gas phase 
is. 



dCO 

dt 



gas 



CO 



gas 



CO 



T„ 
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^off 



(13) 



3.3. CO abundances in model cores 



The combined effects of photodissociation and depletion result in the abundance of CO 
peaking inside the boundary of the core where the visual extinction is just high enough to 
shield the CO, but the density has not yet become high enough for significant depletion by 
freeze-out. This effect is seen in figures |3] and Si Figure [3] shows a spherically-symmetric, 
static-equilibrium model with the properties of a therm ally super-critical core. This core is 
the same as shown in figure 5 of iKeto fc Fieldl (120051 ) but with a factor of 2 higher total 
mass. Figure H] shows a stable, starless core perturbed by an internal oscillation. This 
particular structure of an oscill ating core, motivated by observations of B68, is the same as 
illustrated in figures 6 and 7 of iBroderick et al.l ( 120071 ) and is included here to show how the 
abundance varies in response to an asymmetric density profile. The sharp drop in density 
at the boundary of this thermally sub-critical core is due to the high internal density at the 
boundary where the core is truncated to give it a low mass. In B68, this transition implies 
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a high pressure exterior confining gas, probably hot gas. 



4. Energy Balance and Temperature 
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2005). These models include the effects of molecular line cooling, the ra- 



diative equilibrium of dust in starlight, dust-gas coUis i onal c oupling, and cosmic ray heating. 
We start with the model described in iKeto &: Field! (120051 ) that includes these effects, and 
improve this model by including variations in the cooling rate due to abundance variations 
and by including photoelectric heating at the cloud boundary. 



4.1. Line Cooling with Abundance Variations 



In a model for molecular line cooling that uses cooling coefficients based on a local 
appro ximation such as the large velocity gradient (LVG) radiative transfer model ([Goldsmith 
200ll ). it is straightforward to include the effect of abundance variations due to depletion. To 
do this, we use the cooling coeffi cients for different "depletion factors" , for example as listed 
in the tables in lGoldsmithI (1200 ll ). We assume that the abundance of coolants is given by our 
simple model of the CO abundance. This follows because most of the molecular l ine cooling 
is thr ough transitions of carbon species such as ^^CO, ^^CO, C^^O, C, and CS (iGoldsmith 
200lh . 



At the cloud edge where the CO abundance varies due to dissociation we must also 
include the cooling due to C+ which can be as effective a coolant in photodissoc iation regions 
as CO is in molecular regions. The C^ cooling rate is given by (iTielend (120051 ) , eqn 27), 



Ac+ = 3 X 10"^^n(C+)n(H2)^(l + 0.42(Xi/l x 10"^)) exp(-92/r) ergs cm 



' s-i 



(14) 



where Xi is the ionization fraction, taken to be the abundance relative to H of C"^ 
1.4 X 10-^ C+/Ct 



or 



^total • 



Figures O and E] show the effect of abundance variations on our two classes of cores. 
While a decrease in abundance results in a decrease in the cooling rate, the temperatures 
of the depleted and undepleted cores are about the same. In the thermally super-critical 
cores, the cooling at the core center is dominated by collisional coupling to the dust and the 
cooling at the outer radii is dominated by C^. Only at mid-radii is the cooling primarily 
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through molecular lines. In the thermally sub-critical cores, abundance variations also do 
not change the temperature too much, but for a different reason. The optical depths of 
the coolants in the thermally sub-critical cores are low enough that the gas can always 
find some transitions with moderat e optical depths for ef f ective cooling despite the changes 
in abundance caused by depletion. iPavlyuchenkov et al.l (120061 ) also calculate the effect of 
depletion on the temperature structure of model starless cores. 



4.2. Photoelectric heating 



At the edge of the core, high energy photons can liberate elec trons from dust 



grains 



by the photoelectr ic effect. In cold dark clouds, the heating rate is (IBakes fc Tielenslll994j : 
Young et al.]l2004h . 

(15) 



pe 



10-^\Gpe(r)n(H2) 



ergs cm ^ s ""^ 



0.5 for the conditions in cold dark clouds. The factor Gpe is 



where the efficiency factor e 
the number of high energy photons [hu > 6 eV) normalized by the number in the general 
interstellar radiation field. In computing Gpe{r), the intensity is averaged over frequency 
and direction, f2, 



An 



JydvdVL 



6eV 



Jy exp(— ry(r, uj))dudVL 



(16) 



6eV 



Here Ty{r,uj) is the frequency dependent optical depth from r to the cloud surface r = 
along some particular direction, uj. 

Figures [7| and [S] show the temperature structures of our two classes of cores with pho- 
toelectric heating. Where the gas is thin enough that shielding is ineffective, photoelectric 
heating raises the temperature. Figures [7] and [8] also compare the gas temperatu re with the 
excit a tion temperature of CO computed by our non-LTE radiative transfer code (IKeto et al. 
2004: iKeto fc Field 2005) assuming a total CO abundance relative to H2 of 5.625 x 10~^ 
(jGoldsmithll200ll ). In the thermally super-critical cores, the gas density is very low at the 
boundary (figure [H]), and the CO is not strongly excited by collisions. Thus the excitation 
temperature of CO declines with radius despite the increase in the gas temperature. This is 
consistent with observations of the diffuse interstel lar medium that show H spin ternperat ures 
of about 75 K and CO temperatures below 5 K (jSurgh et al. 2007 : Pineda et al.lboosi ). In 
contrast, because of the high density throughout the sub-critical cores (figure [HD, the CO is 
always coUisionally excited and the CO excitation temperature is everywhere approximately 
the same as the gas temperature. The CO excitation temperature slightly exceeds the local 
gas temperature at locations where the gas temperature is rising rapidly and the CO line 
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has a non-neglible optical depth. This is be cause the hne radiatio n from the nearby warmer 
gas heats the local CO level populations. iPavlyuchenkov et al.l (120061 ) also calculate the 
brightness and ratios of some CO lines in models of starless cores. 

Observationally, the excitation temperature of CO in cores may be determined by com- 
paring the brightness of the (2-1) and (1-0) transitions. Figures [10] and [11] show that the 
ratio of CO (2-1) to CO (1-0) declines at the edge of the thermally super-critical cores as the 
excitation temperature decreases whereas in the thermally sub-critical cores, the ratio rises 
at the edge. This predicted decrea se in the CO excita tion temperature has been observed at 
the edge of TMC-IC (figure 12 of ISchnee et al.l 120071 ). a core whose properties put it in the 
thermally super-critical class. In contrast, at the edge of the sub-critical core B68, the ratio 
of th e CO lines is observed to be constant or possibly slightly rising (figure 5 of iBergin et al. 
20061 ) . Although the decreasing line strengths adversely affect the signal-to-noise just at the 
edge where the comparison is the most diagnostic, the observed ratios are consistent with 
the different behavior of the line ratios expected in each of the two classes of cores. 



5. Implications 

5.1. The rate of cosmic-ray ionization and the opacity of fluffy dust 

Recent observations show promise of measuring the gas temperature with a precision of 
a few degrees. Such precision might allow a determination of the rate of cosmic-ray ionization 
and also the opacity of the dust. In dark clouds and dense cores, observations of molecular 
ions as well as measurements of kinetic temperature limit the possible values of the cosmic- 
ray ionization rate to between about lxlO~^^ and 6xl0~^^ s~^ ( Ivan der Tak fc van Dishoeck 



2000 ; Dalgarno 12006 ). In our models we compare three different rates, 1.3 x 10 s ^ 



( ISpitzer fc Tomaskdll968l ). 3.0 x 10 s ^ and 6.0 x 10 s ^ which we refer to as as low. 



standard, and high cosmic-ray ionization rates. Figure [12] and [13] show the temperature 
structure for these three different rates in our two classes of cores. In the thermally super- 
critical cores, the cosmic rays most strongly affect the gas temperature at mid-radii where the 
energy input from cosmic rays is significant in the thermal balance. In the core center where 
the gas temperature is coUisionally coupled to the dust, the cosmic rays have little effect 
because they carry much less energy than the dust radiation. At outer radii the dominant 
energy source is hot photoelectrically released electrons. The "high" rate pro duces temper 



atures at mid-radi i that are higher than ge nerally indicated by observations (iTafalla et al. 



2004; 


Young et al. 


2004; 


Crapsi et al. 


2007) 



vations at core radii of 0.05 to 0.1 pc are needed to make a definitive statement. In the 
thermally sub-critical cores, because the temperature is nearly uniform except at the edge. 
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and because the gas and dust are not well coupled, the higher cosmic ray rates simply raise 
the gas temperature across the cores. 

The gas temperatu re is also dependent on the dust opacity. From a comparison of CO 



and dust observations, iKruegel fc SiebenmorgenI (119941 ) found that the opacity of dust in 



dark clouds is higher than in the more diffuse interstellar medium. In the cold quiescent 
interiors of dark clouds, interstellar dust grains acquire thick icy mantles that increase their 
opacity. In the high density centers of thermally super-critical cores, t he ice-coated dust may 



coagu late and become fluffy, further increasing the opacity (see also lOssenkopf fc Henning 



(I1994J )). Further evidence for higher dust opacities is provided by observations by Evans et 
al. (2001) that show that a standard dust opacity results in core masses that are beyond the 
gravitational stability limit if the cores were modeled as Bonnor-Ebert spheres. If the dust 
opacity were higher than the standard opacities in the interstellar medium, then the core 
masses would be less, and the cores more stable with lifetimes longer than free-fall times and 
more consistent with the numbers of cores observed. Keto et al. (2004) determined the core 
masses independently of the dust opacity from observations of N2H"'" lines and also found 
indications that the dust mass opacity appropriate for the more diffuse interstellar medium 
is too low for cloud cores. 

Precise measurements of the gas temperature in starless cores could provide yet another 
line of evidence for or against fluffy dust. Figures [T2] and [T3] compare t he gas temperatures in 



our tw o classes of cores calculated with the standard dust opacities of lOssenkopf &: Henning 



( ll994j ) and with the opacities increased by a factor of four. The increased opacities reduce 
the gas temperature everywhere. In the center of the cores, the gas temperature is reduced 
because the radiative input to the dust is decreased. At the outer radii, the increased dust 
opacity provides shielding against the high energy photons responsible for the photoelec- 
tric heating. Although there are complicating factors, the temperatures determined from 
molecular line observations are generally more consistent with the lower gas temperatures 
predicted by the increased dust mass opacities of fluffy dust. This is particularly noticeable 
in the center of t he thermally supe r-critical cores where the observed temperature, about 



6 K, (figure 4 of ICrapsi et al.l 120071 ) is below that obtainable with standard dust opacities 



even with the "low" rate of cosmic ray heating. 

Another way to raise the dust optical depth without changing the dust opacity is to 
raise the gas density. If the density of the core were higher, the dust optical depth would 
also be higher, even for the same mass, and the central temperatures therefore lower. Our 
model assumes a ce ntral density of 1 x 10^ cm~^, a little less than estimated for L1544 by 



Crapsi et al.l (120071 ). In figure [12] we present the results for our lower density to maintain 



consistency throughout our calculations. It is simple enough to compute temperatures for 
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models with higher densities, and the results show that a model with a central density a 
little over 10'' cm~^ is required to bring the central temperature down to 6 K if we use the 
standard opacities and cosmic-ray ionization rate. A higher central density of course implies 
a higher density throughout the core, and this would not be consistent with continuum and 
molecular line observations of L1544. 

Although our models assume a dust opacity that is ind ependent of density, the opa city 



in starless cores may be higher at densities above 10^ cm~^ (jOssenkopf fc Henning]|l994l ). In 
order to put stringent constraints on the cosmic-ray ionization rate, accurate measurements 
of the gas temperature as well as the dust continuum emission at various wavelengths will be 
needed to determine the temperature and dust opacity as a function of radius and density. 



5.2. The dynamical stability of starless cores 



The starless cores are generally successfully modeled as Bonnor-Ebert spheres with the 
addition of some non-therma l energy amounting to a fraction of the total internal energ y 
JSarranco fc GoodmanI llQQsl : ICoodman eTaP llQQsl : iTafalla et all booi Ihada et all boosh . 
Any disagreements with observations seem resolvable by the addition of further physical 
processes to the models. As suggested in our study here, some of these details are significant 
such as non-isothermal gas, variable molecular abundances, and perturbations in velocity and 
density, but so far no observations or theoretical considerations have suggested a revision of 
the basic model that the cores are pressure-supported, self-gravitating clouds. 

However, the Bonnor-Ebert spheres as theoretical objects have requirements that are not 
necessarily met by the starless cores. As truncated solutions of the Lane-Emden equation, 
the BE spheres require an exterior bounding pressure. If the bounding medium is cold 
enough that its own self-gravity is significant, then the combination of the BE sphere and 
the confining medium cannot be dynamically stable. However, if the external medium is hot 
enough that self-gravitational forces are negligible, then the bounding medium itself will not 
collapse onto the sphere and decrease its stability. This seems to be the case for isolated 
cores such as B68. What about cores that exist in the midst of a star forming region filled 
with widespread, low-density, molecular gas? Is the increase in temperature of the molecular 
gas caused by photoelectric heating sufficient to stabilize the surrounding molecular gas? 

To investigate this possibility we calculate the maximum stable central densities for cores 
of different masses with and wi thout photoe lectric heating. This calculation assumes the 
stability criterion of BE spheres (jBonnorlll956l ) applied to non-isothermal cores. As a function 
of mass we compute the maximum or critically stable central density which occurs at the 
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maxim um stable ext ernal pressure as ind i cated by a diagram such as figure 1 of iKeto fc Field 
( 120051 ) or figure 1 of [Lombard &: BertinI (120011 ) . Figure UM shows that photoelectric heating 
improves the stability of larger cores. However, according to figure [T3] the maximum stable 
density is still lower than the central densities suggested by observations. This does not mean 
that observed cores are necessarily unstable. The addition of non-thermal energies such as 
oscillations, turbulen ce, or magnetic f i elds that are not included in this calculation might 



resolve the difference ( GaUi et al. gOOl iKeto fc Fieldll2005l : iKeto et al.ll2006l : iBroderick et al. 



20071 : iHennebelle fc FromangI 120071 ). However, to include this energy in the model we would 



need to know how it affects the equation of state. If the equation of state remains as P ~ p^, 
then the effect on the stability o f the additional non -thermal internal energy is the same as 
increasing the gas temperature (jXeto fc Field 2005). The calculations for figure suggest 
that, consistent with observations, the additional energy need only be a fraction of the 
thermal energy. The analysis suggests that the increased temperature due to photoelectric 
heating is significant in improving the stability of the more massive cores, but not significant 
enough to eliminate the need for some additional internal energy. 



6. Conclusions 

The structure of starless cores has been revisited, taking into account a simple CO 
chemistry and photoelectric heating in the radiative energy balance. The main conclusions 
of this work are: 

1. Observations of the different properties of different starless cores can be understood in the 
context of a simple physical model that shows different structures and behavior depending 
on whether the central density is above or below about 10^ cm~^. This value depends on 
other properties, in particular masses in the range of 1-10 Mq and temperatures around 10 
K, typical values for starless cores. 

2. The low gas temperature (< 7 K) recently observed toward the center of pre-stellar cores 
can only be reproduced if the dust opacity is increased by a factor of a few, suggesting that 
dust grains become fluffy in cold dark cores where the densities are >10^ cm~'^. 

3. The gas temperature of starless cores is quite sensitive to small changes in the cosmic-ray 
ionization rate, (, at volume densities below about 10^ cm~^ (the critical density for gas-dust 
coupling) and at extinction values larger than about 1 mag (where photoelectric heating is 
negligible). Thus, high sensitivity ammonia observations of starless cores can put stringent 
constraints on (. Recent observations toward L1544 suggest ( ~ 1x10^^'' s~^. 

4. Photoelectric heating affects the stability of a core so that it can support a higher central 
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density or a higher external pressure. However the effect is only important for the larger 
starless cores which still require some additional internal energy to be dynamically stable. 
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Table 1. Parameters of models illustrated in the figures 



Figure 
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Photoelectric 


Depletion 


Dust opacity 






ionization rate 


heating 
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low 


off 
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isothermal 


isothermal 
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same as figure 12 
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= 1.3 X 10" 
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& Hennind l|l994 


J); 4x = dust opacities four 
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3r; both : 


= figure compares structures with depleted and undepleted abun- 



dances. 
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Fig. 1. — Variation in the abundances of C, CO, and in the photodissociation region at 
the edge of a cloud as calculated by the simple model described in the text. In this figure, the 
cloud has a constant density of 1000 cm~^ and the interstellar radiation field has a Habing 
flux. Go = 1. 
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Fig. 2. — The relative abundance of CO in the gas phase in equihbrium between depletion 
and desorption as a function of density. The depletion is caused by the loss of CO from the 
gas phase as the molecules freeze onto dust grains (equation [12] in text). The depletion is a 
function of density because the freeze-out is proportional to the dust-gas collision rate and 
inversely proportional to the cosmic-ray induced desorption rate which is itself independent 
of density. 
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Fig. 3. — Structure of a thermally super-critical core. The figure shows the CO abundance, 
mean extinction, and gas density of a model cloud with the parameters of a thermally 
super-critical core. The model cloud is a Bonnor-Ebert sphere similar to that in figure 5 of 



Keto fc Field! (120051 ). The lowest curve is the log of the relative CO abundance. A value 



of zero represents no reduction in CO abundance from photodissociation or depletion. The 
middle curve is the transmission of starlight. A value of 1 represents no extinction and zero 
represents total extinction. The top curve is the log of the number density. Although the 
maximum radius on the plot is 0.6 pc, the core would not be observed to have this extent 
because of the low density and molecular abundance at large radii. Parameters listed in 
table [H 
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Fig. 4. — Structure of a thermally sub-critical core. The figure shows the CO abundance, 
mean extinction, and gas density in the same format as figure [31 Here the model has the 
parameters of a stable, starless, thermally sub-critical core. The asymmetries are perturba- 
ti qns from an internal os cillation. The oscillations are further illustrated in figures 6 and 7 
of iBroderick et al.l (120071 ). Parameters listed in table [TJ 
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Fig. 5. — Comparison of the temperature of thermally super-critical cores, with and without 
variable CO abundances. Although the abundance variations affect the molecular line cooling 
rate, where the abundance is most reduced the cooling is dominated by the dust or C"*" rather 
than the molecular lines. The relative abundance of CO has been multiplied by 10 so that 
a value of 10 means undepleted. Parameters listed in table [H 
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Fig. 6. — Comparison of the temperatures of stable, starless, thermally sub-critical cores, 
with and without variable CO abundances. These cores are more nearly isothermal than the 
thermally super-critical cores (figure [5]) because the gas does not couple with the dust and 
is cooled primarily by line radiation. The relative abundance of CO has been multiplied by 
10 so that a value of 10 means undepleted. While the abundance of CO is quite low across 
the core, there is still a factor of five difference between the center and the edge. Parameters 
listed in table [H 
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Fig. 7. — Comparison of the gas temperature and the excitation temperature of CO (1-0) 
in an unstable, pre-stellar, thermally super-critical core. Because the gas density (figure 3) 
and abundance of CO (figure 5) are so low in the outer core, the excitation temperature of 
CO declines with increasing radius even as the gas temperature increases. In the center of 
the core, the CO lines are optically thick and the excitation temperature is approximately 
the same as the gas temperature. Parameters listed in table [H 
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Fig. 8. — Comparison of the gas temperature and the excitation temperature of CO in a 
stable, starless thermally sub-critical cores. Because of the low mass of this model (1 Mq), 
the core is truncated at a high enough density that the excitation temperature of the CO 
remains close to the gas temperature up to the boundary. Parameters listed in table [TJ 
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Fig. 9. — The density structure of unstable, pre-stellar, thermally super-critical cores, 
compared to the density of stable, starless thermally sub-critical cores. The sub-critical 
cores are truncated at a high density as would be appropriate for low mass cores in a high 
pressure environment. The density of the thermally super-critical cores decreases to a much 
lower value at the edge, consistent with a lower surrounding pressure. Parameters listed in 
table [H 
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Fig. 10. — The line intensity of C^^O(l-O) and C^^0(2-l) (radiation temperature) and their 
ratio (2-l)/(l-0) in unstable, pre-stellar thermally super-critical cores. The temperature 
structure for this model is shown in figure [71 The density structure is shown in figure [91 The 
ratio of the CO lines is a function of the excitation temperature of the CO. In this model 
because of the low densities at the edge of the core, the CO excitation temperature is much 
lower than the gas temperature. Parameters listed in table [H 
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Fig. 11. — The line intensity of ^^CO(l-O) and ^'^C0(2-l) (radiation temperature) and their 
ratio (2-l)/(l-0) in stable, starless thermally sub-critical cores. The temperature structure 
for this model is shown in figure [8l The density structure is shown in figure [91 The ratio of 
the CO lines is a function of the excitation temperature of the CO. In this model because of 
the high densities at the edge of the core, the CO excitation temperature is always close to 
the gas temperature. Parameters listed in table [H 
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Fig. 12. — The gas temperature in thermally super-criticals cores under different conditions. 
The figure shows the gas temperature for a "low" cosmic-ray ionization rate of 1.3 x 10~^^ 
s~\ a "standard" rate of 3.0 x 10^^'' s^^, and a "high" rate of 6.0 x 10"^'' s~^. Because 
cosmic rays heat the gas directly, their rate affects the gas temperature in the intermediate 
region where the cooling is not dominated by dust and the heating is not dominated by 
hot photoelectric electrons. The figure also shows the gas temperature if the dust opacity 
is 4 times the "standard" value. This increase in dust opacity lowers the gas temperature 
everywhere - in the center and mid-radii by lowering the dust temperature, and at outer 
radii by increasing the shielding against high energy photons that generate hot electrons. 
Parameters listed in table [TJ 
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Fig. 13. — Same as figure [12] except that the model here is for stable starless thermally 
sub-critical cores, rather than super-critical. Parameters listed in table [H 
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Fig. 14. — The maximum stable central density of a starless core with (solid line) and 
without (broken line) photoelectric heating at the core boundary. For the mass indicated on 
the abscissa, the log of the maximum density in the center of the core is indicated on the 
ordinate. The figure shows that the inclusion of the photoelectric heating at the core edge 
marginally improves the stability of the more massive cores. The cores here are modeled 
with the only internal energy being thermal. The addition of some non-thermal energy in 
the form of microturbulence or magnetic energy would raise the curves to higher density. 



